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LOW-COST INERTIAL NAVIGATION FOR MODERATE-g MISSIONS 


Shmuel J. Merhav* 
Ames Research Center 


SUMMARY 


This paper describes a low-cost inertial navigation system (INS) concept 
for a broad class of flight missions characterized by moderate accelerations 
and limited attitude variations. Typically, these missions may involve general 
aviation aircraft, helicopters, or remotely piloted vehicles (RPV) . Though 
highly desirable, the use of inertial navigation in these aircraft might be 
precluded by the high cost of INS technology constituting a substantial frac- 
tion of the total aircraft cost. The significance of the moderate accelera- 
tion and limited attitude is reviewed with respect to platform mechanization 
and instrumentation, and a novel hybrid mechanization, partially gimballed and 
partially strapdown, is presented. Implemented by an unbalanced two-axis 
gimbal system, controlled by a two-degree-of-ffeedom gyro, it provides locally 
level two-axis acceleration information, along with pitch and roll measure- 
ments. Heading information is provided by a second gyro mounted in the inner 
gimbal. It is shown that the system error model is equivalent to that of a 
conventional platform with a tilt error determined by the integral of the gyro 
drift rate and an equivalent accelerometer bias proportional to the same drift 
angle. Thus, by calibrating platform drift rate, the accelerometer-type errors 
are also cancelled. Rapid gyrocompassing, implemented with opened gimbal 
control loops, and a strapdown procedure also provides calibration of gyro 
drift rate biases. Being subjected to small angular inputs and a low-g envi- 
ronment, g and g^ dependent errors are negligible. Thus, adequate precision 
can be obtained from moderate cost gyroscopes, resulting in a rms navigation 
error of the order of 1 n. mi, /hr. The dispensing of accelerometers and the 
simplified mechanization and computation imply substantial cost savings and 
improvement of reliability. 


I . INTRODUCTION 


Autonomous navigation, which by definition excludes external navigation 
aids, is essentially based on inertial instrumentation, except for the crudest 
heading and velocity methods using magnetic and airspeed measurements. Origi- 
nally developed for ballistic missiles and military aircraft with all-attitude 
and high-g capability, the technology in a state-of-the-art inertial naviga- 
tion system (INS) of the 1 n. mi. /hr class incurs costs in the $50,000-$100,000 
range. This high cost often precludes its use in a wide scope of applications 
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such as general aviation, helicopters, RPV’s, and ground vehicles. In these 
applications the principal common factor is their low-g environment. It is 
the purpose of this paper to show how these relatively benign dynamic condi- 
tions can be advantageously used in the choice of sensors and IMU mechaniza- 
tion, so as to achieve a substantial reduction in cost and complexity. 

Efforts to reduce cost and size and to improve reliability have, in the 
past decade, been largely oriented to strapdown mechanizations (ref. 1). How- 
®ver, since by far greater demands are placed on the gyroscopes and on 
on-board computational volume and speed, the advantages expected from dispens- 
ing with the gimbal system are largely offset. The greater sensitivity of 
strapdown systems to alignment errors (ref. 2) and to correlated noise compo- 
nents in the accelerometers and gyroscopes (ref. 3) have so far impeded their 
successful competition with gimballed systems in the 1 n. mi. /hr class 
(^efs. 1 and 2). For short-term guidance in tactical missiles or in aided 
navigation with frequent updating (ref. 4), strapdown mechanization is superior 
with respect to miniaturization and reliability. 


Alternative approaches for significant improvement in cost effectiveness 
appear to involve a major departure from classical inertial navigation prin- 
ciples or an entirely new concept in its instrumentation. Examples of such 
departures have recently been reported. By constructing an optimal observer 
(ref. 5) on the (known) parameters of the aircraft, velocities are estimated 
from measurements of roll rate, pitch rate, heading, and altitude. The method 
appears to have good potential for short-term autonomous navigation. However, 
further work is needed with regard to sensitivity to parameter variations, 
trim errors, and wind modeling. An attempt to depart from conventional IMU 
instrumentation by entirely dispensing with gyroscopes has recently been 
reported in reference 6. A single-axis physical pendulum is Schuler tuned by 
measuring the angular acceleration Induced by external specific forces and 
artificially providing the immense moment of inertia of a Schuler pendulum by 
a torque motor, driven by the highly amplified angular acceleration signal. 

Not being gyroscopic, the attractive property of the device is its indiffer- 
ence to the Earth's rotation. However, the excessive gain (~10^) and gain 
stability required in its realization raises doubts as to its performance in 
an actual flight environment. 

An earlier attempt to depart from conventional IMU instrumentation by 
dispensing with accelerometers has been reported by Hector (ref. 7) and Astrom 
and Hector (ref. 8). Instead of a simple pendulum, a gyro-pendulum is Schuler 
tuned, again requiring excessive torquing commands. Additional angular motion 
sensors and torquers are however, required to provide platform isolation from 
vehicular attitude motion. A detailed design and error analysis of this con- 
cept has been reported by Koenke (ref. 9) , verifying that the system obeys the 
basic laws of error propagation in INS. The advantage of dispensing with the 
accelerometers is apparently offset by a large sensitivity to vibration and the 
need for critical adjustments of the large gains. The foregoing examples, 
aimed at reducing the cost of the sensor package, reflect the continuing 
search for low-cost navigation systems, especially for the aforementioned type 
of low-g missions in the 1 n. mi. /hr class for which the technologies of con- 
ventional IMU's, originally developed for military applications, are highly 
overdesigned. 
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In the IMU configuration presented in this paper, the less demanding 
low~g environment is exploited so as to achieve a substantial reduction in the 
number of inertial sensors and in the complexity of mechanization without 
departing from the basic principles of inertial navigation (ref. 5) or resort- 
ing to new instrumentation concepts (ref. 6). It is shown that, with modern 
strapdown gyros, the gimballed and strapdown concepts can be combined into a 
hybrid mechanization, using two gimbals only and dispensing with the acceler- 
ometers. The local-level, two-axis acceleration information is provided by 
the two gimbal torquers, pitch and roll measurements by pick-offs, and the 
heading information by an azimuth gyroscope mounted in the inner gimbal. This 
hybrid mechanization, requiring only the single azimuth Euler transformation, 
which can be executed by a conventional, relatively slow microprocessor, has 
apparently not been considered before. It appears to be superior in cost 
effectiveness to both the gimballed and strapdown mechanizations. 

The feasibility of the method is investigated by modeling the dynamics, 
controls, and error sources in the sensors and system and by statistical 
analysis of the error propagation for flight durations of 1 hr. If experimen- 
tally validated in the future, the proposed solution might lead to a substan- 
tial extension of autonomous navigation to the types of aircraft in which low 
cost is of paramount importance. 


II. SYSTEM CONSIDERATIONS 


In this section, INS mechanizations are reviewed in the light of the 
assumptions and requirements characterizing moderate-g missions. The conclu- 
sions point to the hybrid mechanization presented in this paper. 


A. Assumptions 

The class of mission profiles, for which low-cost autonomous navigation 
is considered in this report, is characterized as follows: 

1. Atmospheric flight at moderate altitudes and velocities typical of 
general aviation, helicopters, RPV’s and surface vehicles. 

2. The specific forces acting on the vehicle do not exceed the order of 

1 g. 

3. Unlimited attitude except in pitch, which is limited to ^±80°, as 
typical in vehicles controlled by vertical gyroscopes. 

4. Yaw rate is assumed not to exceed values on the order of 10° /sec. 

5. Navigation is implemented in conventional spherical coordinates, dis- 
regarding Earth oblateness. 

6. The magnitude and direction of the gravity vector are known constants 
throughout the flight mission. 
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7. Duration of a typical autonomous flight phase does not exceed ~1 hr. 

8. Altitude information is derived from conventional altimeters (baro- 
metric or other) . 

9. The navigation error defined in the local-level plane during the 
autonomous flight phase should not exceed the order of In. mi. 


B. Review of Mechanization Problems 


1. Navigation Equations 

The relation between acceleration R, specific force and gravitation 
G in a reference frame rotating at an angular rate oj with respect to iner- 
tial space is 


F + G= R+ cjxR+2a)xR+u)X (ojxR) 


( 1 ) 


Expressed in geographic coordinates of a perfectly spherical Earth, equa- 
tion (1) takes the form of the nonlinear time-varying differential equation 

+ ^ + r(^, X, R) (2) 


where ^ is the net gravity vector given by 

"O" 


S. = 


0 

L8J 


= G - n X (Q X R) 
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(3) 


and is a known vector valued function of Vq, and the scalar parameters 

fis, X, R. Explicitly, equation (2), with equation (3), takes the form 
(ref. 3) : 
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in which X is determined from the integration of X = Vjj/R. Vjj, Vg, Vy, Fj^, 
Fe> Fv defined in the local level, north, east, and vertical coordinate 

frame. Any mechanization of INS requires: 
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1. Determination of %» ^E’ ^V 

2. Integration of equation (4). 


In view of assumption 8, Vy can be derived from the independent measurement 
of altitude h, so that equation (4) reduces to the two-dimensional navigation 
equation, in which Vy is now a time-varying parameter 





2fig sin X + 
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sin X + tan X]V^, 4* cos X + 


V V 
E V 


(5) 


Measures of Fg, are normally obtained from onboard accelerometers 

mounted in a coordinate frame (X^, • The measured specific force vec- 

tor = col[Fxa> ^ya> Fzal tnust be transformed to yield F^ and Fg. The 
three fundamental aspects of inertial navigation emerging from equation (5) 
are 


- Instrumentation by which the specific force vector _F is measured in 
an airborne coordinate frame defining F^ . 

- Transformation — by which these measurements are transformed into the 
desired local-level geographical coordinate system, = Zc/aZa* 

- Navigational computation — by which equation (5) is integrated to yield 
Vn, Vg and, by further integration, vehicle position Xq, Yq. 

2. Instrumentation 

Specific Force Measurement- In view of assumption 2, the required measure- 
ment range is 1 g and by assumption 9, the allowable null point uncertainty 
should not exceed 100 yg (appendix C) . Thus, the dynamic range is of the order 
of 10^ instead of 10^ in military applications, and 10^ or better in long-range 
missions such as transatlantic flight or submarines. This reduction in dynamic 
range is a key feature in the realization of the low cost of the IMU described 
in the following section. 

Platform Orientation Error- Assuming that the orientation errors in gim- 
balled and strapdown mechanizations are equal, the error of the resulting 
specific force measurement is determined as follows: For unaccelerated motion 

along the vertical axis Vy = 0 in equation (4), the specific force vector 
acting on the vehicle is ^ = co1[Fn, Fg, -g] . An orientation error vector 
^ = co1[4j^, causes the error AF^g col[AFig, AFg, AFy] . 
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By equation (5), only AFfj, AFg are relevant. Thus, equation 
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The first term in equation (7) is the d is tance ~r elated navigation error, pro 
portional to the azimuth error The second term is the time- related navi- 

gation error, due to the platform tilt with respect to the local vertical. 

Required Performance of Gyroscopes- In view of assumptions 1 and 7, for 
an aircraft with a ground speed V = 300 m/sec, and a flight duration of 1 hr, 
the range covered is 1080 km* For a 1 n. mi* error, 4 'a ~ can be tolerated. 

= 0*1® /hr would result in an error of 0.5 n. mi* only. and i|^£, however, 

must not exceed 0.01® /hr as indicated in appendix A. In the gimballed ^mecha- 
nization, the tolerable error in is thus '^lO times the error in and 

In the strapdown mechanization, however, all three angular rate measure- 
ments interact to yield the required direction cosine matrix. Thus, all 
angular rate measurements must be in the 0.01® /hr class. 

Furthermore, the rate input range of the gyros is directly dictated by 
the attitude rates of the aircraft. As a result, the dynamic range 
of a strapdown gyro for missions characterized in assumptions 1 through 9 is 
from 0.01®/hr to 100®/sec, i.e., a dynamic range of 36x10 , as compared to 
gimballed platforms in which the required dynamic range is from 0.01® /hr to 
10®/hr (i.e., a dynamic range of 10^). The large input rates in strapdown^ 
systems, based on spinning wheel type gyros, incur intensive torquer activity, 
which causes greater temperature fluctuations and reduced accuracy. Also, the 
excessive demands on torquer scale-factor accuracy essentially exclude 
moderate— cost spinning-wheel gyros from strapdown systems, even for the rela- 
tively benign mission requirements defined in assumptions 1 through 4. 

Alignment and Calibration- A comparison of the effect of gyro and acceler- 
ometer alignment error on navigation accuracy for gimballed and strapdown sys- 
tems is given in reference 2. It is indicated that the analytical leveling in 
the strapdown mechanization is valid only for the airframe orientation on the 
ground at which the leveling was performed. In subsequent turns during flight, 
the computer— derived axis system tilts with respect to the instrument axis sys- 
tem. This effect does not exist in gimballed systems. Gyro axis misalignment 
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also causes by far a larger error in strapdown systems than In gimballed sys- 
tems. These errors are flight-profile-dependent and cannot easily be detected 
on the ground. A special fixture for tilting the entire strapdown system 
through precise angles is required for calibration and to cancel flight- 
profile-dependent calibration errors. 

Sensitivity to Dynamical Environment- The mechanization of XjS/h 
determined by gyro measurements and by the accelerometers. Mechanical 

vibration of the instrument package may exite correlated signals in and 

Fb and cause a bias in £q. This susceptabillty does not actually exist in a 
gimballed system. 

3. Transformation 

In the gimballed mechanization with its coordinate frame (Xp, Yp, Zp) 
perfectly aligned with (Xq, Yq, Zq) , the transformation of 

I^p = col[Fxp, Fyp, Fzp] to is = Tp/p = I and ideally _Fp = Fg. In 

view of assumption 3, the mechanization normally requires at least three 
gimbals . 

In the strapdown mechanization, the instrumentation axis frame being the 

aircraft frame body axis (Xb, Yb, Zb), lb = col[Fxb, Fyb> ^zb^* Ic/a = Ic/b 
and 


4 


^/b^b 


( 8 ) 


where 


Ic/i - [♦l.lelyWl, (9) 

is the nine-component direction cosine transformation matrix. It is self- 
evident, that due to the all-attitude capability (assumption 3), all three com- 
ponents Fxb, Fyb, Fzb must be measured to determine Fjj, Fg. The mechaniza- 
tion providing Ic/b obtained from the following differential equation 
(ref. 3): 

^/b ~^/b^b -^^G^/b 


The skew symmetric 3x3 matrix consists of elements determined by 

the three body axis angular rate measurements p, q, and r; the skew symmetric 
matrix consists of the three angular rate elements of the rotating coor- 

dinate geographic frame, at the location of the vehicle. It follows that 
since the angular rate measurements p, q, and r and all three components of 
Xb are required, assumptions 1 through 9 do not permit savings or simplifica- 
tions in the strapdown mechanization. 

4. Computation 

In the gimballed mechanization, the computation required is essentially 
the integration of equation (5). Since the second term consists of compara- 
tively slowly varying parameters and variables and in view of assumptions 1, 
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2, and 4, they will be three orders of magnitude smaller than the first term, 
the computation of the second term can be performed at an iteration rate in the 
range of 1 to 10 Hz. The first term is normally integrated at a 40-Hz itera- 
tion rate. In the strapdown mechanization, in addition to integrating equa- 
tion (5), equation (10) must be integrated and equation (8) must be computed 
at a normal iteration rate of 40 Hz. Thus, assumptions 1 through 9 do not 
lead to compromises with regard to computational volume and speed. A compari- 
son of a completely coded navigation program for a strapdown system with a 
program for a gimballed INS shows that the latter uses 30% of the real-time 
computer capacity and 30% less memory storage than does the former (ref. 2). 


C. Summary 

The foregoing comparison of gimballed and strapdown mechanizations, based 
on low-cost spinning-wheel gyroscopes, clearly indicates the superiority of 
the former for the class of missions considered in this report. The main con- 
clusions, as they relate to assumptions 1 through 9, are summarized in table I. 


D. Rationale for a Hybrid Mechanization 

The foregoing system considerations and conclusions explain the limita- 
tions of a cost-effective wide spread implementation of a strapdown or a gim- 
balled INS for the class of missions considered here. However, by combining 
the most essential functions of the gimballed concept (isolation of the level- 
ing gyro) with the least critical function of the strapdown concept (azimuth 
Euler transformation) into a hybrid mechanization, major savings in cost and 
complexity can be achieved. This hybridization is made feasible by the limited 
attitude requirement, the low-g environment, and the availability of strapdown 
type gyros (large angular rate inputs). 

The hybrid mechanization, as described in the next chapter, consists of a 
two-gimbal system and two gyros only. Acceleration information is provided by 
the gimbal torquers and heading information from the azimuth gyro which, being 
mounted in the inner gimbal, is exposed only to the relatively low yaw-rate 
inputs and still meet the less stringent requirements for ifi as Indicated in 
section II B. 

The dispensing with the accelerometers and one gimbal, while essentially 
retaining the computational simplicity of the gimballed mechanization, consti- 
tutes a major saving in cost. 


III. DESCRIPTION OF THE GYROSCOPIC NAVIGATION SYSTEM (INS) 


A complete qualitative description of the hybrid mechanization and INS 
integration is given in this section, along with definitions of axis systems 
and symbols used in the analyses of the next two sections. 
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A. The Gyroscopic Inertial Measurement Unit (GIMU) 


The device, schematically shown in figure 1, consists of the following 
elements: a leveling gyroscope LEG; an azimuth gyroscope AZG; and a two-axis 

gimbal system g^, g^, with its torquers T^j Ty, and pick-offs P<j, and P0. 

LEG and AZG are assumed to be identical and of the strapdown TDF type. The 
gimbal assembly is essentially of the type used in conventional vertical 
gyroscopes with relatively powerful torquers T^j Ty, low friction synchro-type 
pick-offs P|j, and Pe, and low friction bearings b - b*. The outer gimbal go 
is hinged to the aircraft frame AF by means of the bearings b - b’ along the 
longitudinal body axis Xjj. LEG and AZG are mounted in the inner gimbal g£ 
and together with an optional bob weight W, impart an intentional mass unbal- 
ance to the two-axis gimbal system with its center of mass C.G., Z cm below 
the suspension point, that is, the intersection of axes X i, Ygi. The gimbal 
axis Zgj_ tends to align itself with so that n and | tend to zero. 

The specific force vector JF is defined here by Fy£, in the local- 

level orthogonal coordinate frame X^, Y^, Z^. is defined along the pro- 
jection of Xb on OX^Yj^. It is assumed here that = ~g. (This assump- 
tion can be removed, as shown in section IV.) and Fy£ excite n and 

which are sensed by the LEG. The corresponding outputs and ampli- 

fied and integrated in the networks Gy, G^ (detailed in fig. 4). The outputs 
of Gy and G^^ are applied to the gimbal torquers T^ and Ty, counteracting 
the reaction torques induced by F^^ and Fy^. With integral control in Gy 
and Gx. n and ? essentially are zero in the steady state. The outputs of 
Gy and are proportional to and Fy£, and are denoted by F^p, Fyp, 

respectively. does not exert a torque and is not measured. The dynamic 

range of F^p, Fyp is determined by the ratio of the maximum torques available 
from T^ and Ty to the friction torques Tf in the gimbal bearings. Realis- 
tic values of these torques show that a dynamic range of ^10^ can be obtained. 

The device performs as a two-axis, horizontally stabilized accelerometer 
independent of vehicle roll and pitch 0 but dependent on ijj. Assuming 
that the ^ before takeoff is known, any subsequent yaw rate sensed 

by AZG is integrated by In\|;» ^ is continuously determined. Therefore, 

this device constitutes a complete IMU and will further on be referred to as 
the Gyroscopic Inertial Measurement Unit (GIMU). It does not require accel- 
erometers and has only two gimbals. This simplified configuration constitutes 
significant savings in complexity and cost in comparison to both the gimballed 
and strapdown mechanizations listed in table^I. Since, in view of section II, 
the accuracy in is '^O.l^/hr and since ^niax “ 10° /sec, the dynamic range 

required from the AZG is only 3.6^10^, which is '^lOO times smaller than in 
typical strapdown applications. A low-cost gyro can meet this requirement. 

From the criteria of section II, the LEG must have a null point stability of 
0.01°/hr. However, this gyro essentially operates in a gimballed mode and its 
required dynamic range is only 10^. With proper temperature control and 
g-compensation, the required null stability can apparently be achieved even 
with a low-cost gyro (appendix E) . This problem is treated in section VI. 

The pitch (0) and roll ((})) measurements, obtained from Pq and P^, are 
insensitive to acceleration and can be used for flight control and display 
purposes in a conventional manner. 
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B. Integration of the GIMU into a Complete INS 

Ip = col[Fxp, Fyp] and rl> provided by the GIMU, the transformation 
defxned by equation (8) reduces to 


= [T]^F 
z— p 

and is specifically implemented in the navigation computer by: 

Ff 


( 11 ) 




N 




L eJ 



cos Ip -sin :p 

sin ip cos ip 


xp 


iJypJ 


( 12 ) 


With Fn, Fg, equation (5) can be integrated to yield and Vg and by 
further integration, the position Xq and Yq. The following gyro precession 
commands, determined in the navigation computer and fed to the gyro torquers 
constitute the complete INS: ’ 

^chuler tuning: derived from Vjj, Vg in terms of body 

axes, fed to T^. 

Earth rate co mpensation : 'I'c^ derived from and X fed 

Tg. T^- 

3. Compensation of g-dependent gyro drift : Derived from the known 

g-sensitivity coefficients of the gyros and from the actual values of F 
Fyp fed to T^, T^, T^. • 

Compensatio n of drift rate bias : derived from a cali- 

bration procedure and fed to Tj^, T^, T^. 

The foregoing compensations, as treated in detail in section V, indicate 
that the complete INS based on the GIMU is feasible. The computations 
involved are straightforward and can be implemented by a conventional 
microprocessor . 


IV. ANALYSIS OF SYSTEM DYNAMICS AND CONTROL 


In this section, the mathematical model of the gyro— pendulum assembly 
shown in figure 1 is developed. The resulting equations of motion are then 
used in the analysis of the closed loop GIMU. 


A. Open-Loop Gyro-Pendulum Dynamics 

The LEG and AZG, assumed to be of TDF DTR type, are described in detail 
in reference 10. Their open-loop transfer functions are equivalent to those 
of a classical two-axis symmetrical free rotor gyro 
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(13) 


0x(S) = -<|)x(S) 


+ 


SMx(S) 


I’^S 



0y(S) = -((.y(S) 


M^(S) + SMy(S) 


I’='S 


+ 


W] 


(14) 


where (x, y, z) is the gyro-casing reference frame, (j)^, (j)^, the angular inputs, 
resolved along x and y; 0x and 0y, are the angular outputs; that is, deflec- 
tions of the rotor with respect to the gyro case; and My are the torques 
applied to the rotor, resolved along the gyro case coordinate frame; H is the 
angular momentum of the rotor; is the moment of inertia of the rotor about 
its X and y axes. The loop closure of the gyro is implemented by setting 

My(S) = A(S)0^(S) (15) 


and 


K^(S) = -A(S)0y(S) 


(16) 


By proper design of the torquing amplifier A(S), the closed-loop trans- 
fer functions are essentially 




03^(S) = -(H/A(S))(J)^(S) 

(17) 

and 






0y(S) = -(H/A(S))$y(S) 

(18) 

and 

from equations (15) 

and (16) the torques are 




My = 

(19) 



I 

II 

(20) 

For 

the LEG mounted in 

the inner gimbal in accordance with figure 1, My 

= -Hi 

and 

Mx = -Hf|. The AZG 

may be mounted with its spin axis collnear with 

either 


the Xgi or the Yg^ axis. Assuming the former, its reaction torque 

= -Hn is exerted on the complete airframe and has negligible effect. The 
reaction torque due to the airframe yaw rate ip, is^ My = Hi{j is disregarded in 
this section in which ip is assumed to be zero. HiJ; is treated as a torque 
disturbance in section VI. 
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The development of the equations of motion of the gimbal system is based 
on reference 11, For small values of rj, the equations reduce to 


where : 


it + Hn = -Wil5 

jfi - h| = -W5,n 


1 = 1^ + \ 

X ^ 


( 21 ) 

( 22 ) 


J=Iy 


H = n 


r r r 
I = I = I 

-Ly j-x J- 


X ■ ‘X 

gi 


(23) 


With specific force components and Fy£ defined along and Y^, equa- 

tions (21) and (22) are modified to include the corresponding inertia torques 
(see fig. 2 for sign conventions) 



(24) 


(25) 


Dividing equations (24) and (25) by I and J, respectively, we have 

C + an = -y(g5 - Fy^) (26) 

n - = ”V(gn + F^^) (27) 

where a = H/I, b = H/J, y = m£/l, v = m£/J. Laplace transforming equa- 
tions (26) and (27) and rearranging results in 

[(S^ + yg)(s2 + vg) + abS^lC = (S^ + vg)yF^^ - aSvF^^ (28) 

[(S^ + yg)(s2 + vg) + abS^ln = -(s^ + 

The characteristic equation is the quartic 

s'* + tg(y + v) + ab]S^ + yvg2 = 0 (30) 

Its solution is 


c2 [g(y + v) + ab] -i /[g(y + v) + ab]^ ? 

'’ 1,2 2 -y 4 


(31) 
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It is readily verified that the discriminant in equation (31) is nonnegative. 
Therefore the two roots are 


and 


= +j |.[g(p + v) + ab ] j /[g(y + V) + ab ]^ 

= I [g(y + V ) + ab ] ^| /[g(y + V ) 


1/2 


= ±j03i (32) 


= ±j03, (33) 


Equations (32) and (33) indicate an undamped conical oscillation at the two 
frequencies and 1 ^ 2 - As an example, the following numerical values are 
assumed 


H - 8 gem sec ; I = 1 gem sec^ ; J = 0.5 gem sec^ ; m£ = 0.5 g sec^ 

Thus, a = 8 sec“^; b = 16 sec“^; y “ 0.5 cm“^; v = 1 cm“^; and 
g = 981 cm/sec^. From equations (31) and (32) 

o)j = 20.1 radsec“^ 

0)2 = 34.6 radsec“^ . 

Expressing equations (28) and (29) in terms of and co^ in 

response to step inputs of and F„£, the solutions for ?(S) and n(S) are 


?(S) = 


y(S^ + vg) 


"y£’ 

F. 


av 


n(s) = - 


S(S^ + Wj^)(S^ + 0)2^) (S^ + 0)^2) (S^ + 0)2^) 

by 


'x£ 


v(S^ + yg) 


S(S2 + cJj2)(s 2 + jj^2) 

The solution in time domain is 


■ x£ 


(S2 + 0)^2) (s2 + 0)2^) 


(34) 


(35) 


C(t) = (to, 2 


N vg - to. 


2 vg - t02^ 

- + 

2 

0)2 


\ ^ “ 1 ' 

) - M~:;r 


cos (A)^t 


vg “ 03. 


) /sin 03^t sin W 2 ^\ 


(36) 
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and 



In the n phase plane, equations (36) and (37) represent a two-mode 
elliptical trajectory at frequencies coj and w, crossing the origin and having 
the mean values 5 and n 



(38) 


(39) 


Equations (38) and (39) show that the sensitivities of ? and n to Fyjj^ and 
equal. A qualitative representation of the solution, equa- 
tions (36) and (37), is shown in figure 3. 


In the foregoing analysis of the pendulum dynamics, frictional damping 
has been disregarded. Its effect is qualitatively indicated by the dotted 
trajectory in figure 3. 


In the absence of F^jj, and Fy/j,, a similar elliptical coning motion with 
its center at the origin, 0, would occur. It is determined by setting to zero 
the right-hand side of equations (28) and (29) and determining the initial 
conditions 5(0), 5(0), . . ., and n(0), n(0) .... 


B. Closed-Loop Analysis 


,The loop closure of the pendulum is schematically described in figure 1. 
If 5jj and Hd denote the unknown drift rates of the LEG, then the outputs are 


m 


m 



(40) 


Considering initial tilt angles 5i and ni» the closed— loop equations, by 
equations (24) and (25) are 


15 + Hf, - + 6j) - Fy^] - 


(41) 
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(42) 


Jn - h| = -mil[g(Ti + - Gybyf|j„ 


bx and by are the torque coefficients of the gimbal torquers T^, Ty. In 
order to implement integral control, Gx, Gy are: ^ 


= >=., + s . 

fit + >^3 //<"!' 'it 

(43) 


(’dt + Ky^JJdf dt 

(44) 


Substituting equations (40), (43), and (44) into equations (41) and (42) 
dividing through by I, J, respectively, Laplace transforming and rearranging, 
the solutions of 5(S) and n(S) are 


where 


C(S) 


X(S)Ay(S) - S^aYCS) 
A^(S)Ay(S) + abS^* 


Y(S)A^(S) + s2bX(S) 

Ajj.(S)Ay(S) + abS't 


(45) 


(46) 


X(S) = S(iiFy^ - pgq) + (Ax^s2 + Ax^S + A^^)l^/S (47) 

Y(S) = S(-vF - Vgn.) + (Ay s2 + Ay S + Ay )n,/S (48) 


Ax(S) = S3 + A^^s2 + (A^^ + 

Ay(S) = S3 + Ay^s2 + (Ay^ + 

i 1=1. 

Ay. S KyA^.! * 1 ° 1’ 

Since Ay^ are large compared 

sions (45) and (46) simplify to 


yg)S + A^ 

(49) 

vg)S + Ay^ 

(50) 

2, 3 

(51) 

2, 3 

(52) 


with a and b, the sixth order expres- 


Sy(F - gc.) + s2 + A^ s + A )L/S 
^(S) = X(S) ^ ^ ^3 d 

$3 + + vg)S + A^^ 


(53) 


15 


(54) 


n(S) 


^ gy + ys -h Ay^)n^/s 


Ay(S) 


s^ + 


AyjS^ + (Ay^ + vg)S + Ay^ 


• , 

Since and are assumed to be slow processes, equations (53) and (54) 
reduce to the following for S 0: 


?(t) ^ /ij dt 

n(t) - J 

That is, the gimbal system tracks the LEG gyro tilt drift angles, 
expected. 

The outputs from and Gy, that is, the torquing signals of T^ 
are ^ ^ 

“x = Gxim = 


(55) 

(56) 


as 


and Ty 

(57) 

(58) 


Differentiating ^(S) and n(S) in equations (53) and (54), substituting the 
result into equations (57) and (58), rearranging and substituting Kx. and 
%1 from equations (51) and (52) , the results are ^ 


S2 + A„ S + A^ 

I ^1 ^2 ^3 

‘^x S3 + Ax S2 + (A^ + vg)S + A ^ ^ ^ 

1 2 -^3 


niy 


s3 + Ay^s2 + (Ay^ + vg)S + Ay^ 


[-vFy^ - vgn^ - vgf|^/S - Sq^] (60) 


With knowledge of the nominal values In, Jn> and bxn» byn, the measured spe- 
cific forces Fyn and F^n for S 0 are defined by Fyn = m^by /!„ and 
Fxp = myby^/Jn. Specifically. x Xn n 


^xp = + sn^/s) 


(62) 


The coefficients 3y - (l/bx)(bx /in) 6x ~ ('7/by)(by /J^) ideally are 

equal to unity. 
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The closed-loop system, thus, is equivalent to two orthogonal accelerom- 
eters with scale factor errors 1 - Bx and 1 - By, mounted on a platform 
with the^ initial tilt errors and and the corresponding tilt drift 
errors and f|(j/S. It is a unique property of this platform that acceler- 

ometer bias or alignment errors are nonexistent, since the instrumentation 
error stems from the single source f|^. 

Typical values for and A^^ , yielding closed-loop poles at 

“l,2 “ “3 “ = 0.5x981 = 490 are: = 100, 

^x^ “ 2,810, and Aj^^ = 36,000. Similarly, for the same closed-loop poles, 

for vg = 1x981, Ay = 100, Ay = 2,320, A„ = 36,000. 

•^1 ^2 ^3 

These values provide a bandwidth of ~5 Hz, which is adequate for inertial 
navigation. 


A block diagram describing the closed-loop system is shown in figure 4. 
In it, the transfer functions in the gyro-pendulum assembly, as derived from 
equations (28) to (35), are given by 


k(S) 

n(s) 

The switches 


_ )jiS(S^ + vg) p avS^ 

(S2 + a)j2)(g2 + (j^2) (g2 + ^^^2) (g2 + ^^^2) 

vs(S^ + pg) ^ ^ byS^ 

(S^ + 03j^)(S^ + ( 1 ) 2 ^) (S^ + a)j^)(S^ + u^^) 

Sx and Sy in the torquing loops are explained in 




jy, (64) 


section V. 


V. SCHULER TUNING, GYRO COMP AS SING, AND CALIBRATION 


A, Implementation of Schuler Tuning 


^ Schuler tuning is assured by applying the required precession commands 
kc f)c to^the torquers T^ and T^ of the LEG. These commands are 

obtained from 5 e “ respectively (see^fig. 5a), where Vg 

and are determined by the solution of equation (5). and must be 

implemented in the X^, Y^, coordinate frame, shown on the horizontal pro- 

jection in figure 5(b). The corresponding transformed components are 


and by setting 


• s 

^c 



cos 
sin ^ 


-sin 
cos ^ 


-Vj,/R 


(65) 



the required precession is implemented . 
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For the short-time and moderate-speed missions considered here as defined 
in section II, the analysis of the Schuler tuning can be based on the 
approximation that the second term on the right-hand side of equation (5) is 
negligible (ref. 12), that is. 


V = F 
N N 


V ^ F . 
E E ) 


(65) 


Thus, small amplitude components at the sidereal frequency fig are disre- 
garded. Vn and Vg are obtained by integrating equation (12). In accordance 
with equation (7) , the tilt angles and are substituted for and 
in the X^, Y^, frame, so that the measured specific forces F~ and F„„ 
in equation (12) are ^ 

^xp “ ^xZ (67) 

and 

Pyp “ Fyji - g?° (68) 

where, in accordance with equations (61) and (62), 

n® = + fi^/S (69) 

and 

r = ?i + l^/s (70) 


It is assumed here that 3^ “ By ^ 


From equations (65), (67), and (68) the required precession rates in the 
X£, Y^, Z£ axis frame are 


The analysis of equation (71) must be performed in geographical coordinates, 
namely 


1 

R 


cos \p -sin 

sin if; cos 


in If; ^ c 
s if; J s 


cos if; -sin if; 

sin if; cos if; 


•xp 


'ypl 


dt 


(71) 


■■"-I 1 fP 


cos if; -sin if; 

sin if; cos if; 


xp 


dt 


■ypl 


(72) 


With torquing commands applied to and T^, Fyp and Fj^p are determined with 
respect to the rotating local vertical, that is. 


18 



and 




( 73 ) 


^yp = ^yi - 85° - 8(? - 


(74) 


Substituting equations (73) and (74) into equation (72) and 
Euler transformation, equation (65) yields fiw = -Vvr/R and 
we have ^ 


performing the 
5 e = Vg/R, 


R + g J* n® dt + g J dt + g JJ dt' dtj 

and 

■ R ‘•t - S / C° dt - g / fg ® 'it) 

expression the small angles ? and n (dropping Indices) are 
treated as vectors, so that ^ 


(75) 


(76) 


rijj = n cos 1 |) + 5 sin :p 
“Cg = -n sin \p + ^ cos ip 


(77) 


Laplace transforming and rearranging equations 
the approximation of equation (66), 


(75) and (76) gives, in view of 


Sr, = - ^ _ -S. „ 

'n r rs ’^n 


V 

O _ -8. g N 

RS " r2 s2 


V. 


— -S- r° _ r ^ -g_ E 

E R RS RS ^E j^2 g2 

The solutions for nj^ and Cg are, with g/R = m 

s 

2 


n„(S) - - 


03^ 


+ 03, 


-n° ^ 

2 N RS 


(78) 

(79) 


(80) 


5i.(S) = - 


0)e 


s2 + 0) 


— ^ 

2 RS 


(81) 


Equations (80) and (81) reveal the familiar sinusoidal Schuler frequency 
oscillation due to hn and and the rotation resulting from vehicular 
velocities and Vg. n§ and Cg are functions of \p, since they originate 
rom errors in the X^, Yj^, axes. Thus, changes in ip during flight 
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modulate the random process and f|j, so that the error propagation also 
depends on the actual flightpath. 


The actual specific force components, Fjj, Fg, in geographical coordinates, 
are obtained by transforming equations (73) and (74) into geographical coordi- 
nates, using equation (66) and substituting equations (80) and (81) 


+ g(- 


N 


S2 + „ 2 "n RS Es) - \ + 8 




— n ° 
2 N 


"e = ''e - 


\ s2 + a.32 ^ Rs) E 8 ^ 


2 "E 


(82) 


(83) 


The description of the Schuler tuning implementation is incorporated in 
figure 6. 


B. Gyrocompassing 


The mass unbalance of the gyro pendulum is advantageously used to achieve 
rapid gyrocompassing. On the ground, with the switches S^, Sy (as shown in 
fig. 4) open, the pendulum theoretically aligns Itself with the local vertical, 
as predicted by equations (36) and (37) with F^£ = Fy^ = 0. The deviation 

from true vertical is determined by two factors: 1) the friction torques of 

the gimbal bearings, and 2) the reaction torques of the gyros due to J2g. 

!• For a typical bearing friction torque of T£ = 0.05 g cm and a mass 

unbalance of W£ = 200 g cm, the deviation or Aq is 

^f 0 05 

" wl " Too" " 

2. The reaction torques exerted on the gimbal system by the LEG and AZG, 
due to the Earth rate components JJg cos X and fig sin X, respectively, are 
nearly equal at moderate latitudes (X ~ 45°) and amount to ~10°/hr. For an 

angular momentum of the gyros where H = 8 g cm sec, the total reaction torque 
is 


T = 2 X 8 X 10 

r 3600 X 57.3 


8 X 10 ^ g cm 


This is negligible in comparison with the friction torque and can be disre- 
garded. The precision requirement in azimuth measurement is determined as 
follows: At a heading angle ip, the Earth rate components in the Xjj,, Yj,, Zjj, 

frame are 
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( 85 ) 


“ - 


- 


» m 



^sx£ 


COS Ip sin 4; 0 


COS X 


cos rp Qg cos X 

^sy£ 

= 

-sin ip cos ^ 0 


0 

= 

-sin \p fig cos X 

^szJl 


0 0 1 


sin X 


fig sin X 

- 


- 






As a result of the orientation error vector col [A?, An, At|;], the errors in 
sensing the Earth rate components are 


• 




- 


- 



- 


- 


AC 


0 

Alp 

-An 


COS ^ 


COS 

X 


sin \p 

fig cos X - Anfig sin X 

An 

= 

-A4; 

0 

AC 


-sin ip 

^s 

cos 

X 

= 

-Alp cos ip 

fig cos X + ACJ^g sin X 

Ai|i 


An 

-AC 

0 


^s 

sin 

X 



An cos ip 

fig cos X + AC sin cos X 

“ - 


- 


- 


- 



- 


- 
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In view of equation (84), and since fig sin X and f2g cos X are of the order 
of 10°/hr, the elements multiplying An or A? are of the order of 0.0025°/hr 
(or less) and can be neglected. The remaining terms are essentially 


1 

1 ^ 


-Ai|» 

sin ip 

fig 

cos X 

An 

= 

-Aif) 

cos Ip 

^s 

cos X 

A^ 



0 



The largest tolerable errors a|. An in the 1 n. mi. /hr class 
Thus, from equation (87) 


(87) 


is O.OlVhr. 


Aijj = 


0.01°/hr 
10° /hr 


0.001 rad = 0.057° 


Therefore, gyrocompassing must be accomplished with a precision of at least 
0.057°. 


The longitudinal aircraft axis on the ground deviates from the geographi- 
cal north by the unknown angle \p. The measured rate outputs of the LEG are 
then 


— “ 


— 

— 


— — 


— - 


- - 



5m 


cos Ip 

sin Ip 


cos X 




cos Ip 



• 

m 


-sin ip 

cos Ip 


0 

— 

y. 

— 

-sin Ip 

Qq cosX - 

"d 

^ J 


- 

- 


— — 




. — 




( 88 ) 


• , 

The values for ?jjj and Pm are fed to a resolver with ^ a variable input angle 
as indicated in figure 6. The resolved outputs | and f| are 

K. iv 
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cos(iJ; - lij.) - cos cos X) + sin il; . (fi ,/J^„ cos X) 

1 X Cl s ± a ^ 


-sin(ijj - \p.) “ cos cos X) - sin cos X) 

1 1 d ^ 1 d ^ 


Qg cos X 


(89) 


}\). is varied until f) = 0, so 
1 K 


sin(i{» - ij;.) = \p - [p, = All) = -sin ip.(k,/Q„ cos X) - cos cos X) 

1 1 id° xd^ 

(90) 

Denoting 5d/^s X = a; X = b; and a/ (a^ + = cos a; 

h/(a^ + b^)^/^ = sin a we have 

sin(i|; - ip^) = -(sin cos a + cos \{;^ sin a) (a^ + b^)^/^ 

= -sin(ijj^ + a) (a^ + b^)^/^ (91) 


Example ; 


a = 0,03°/hr/10°/hr = 0.003 
b = 0.02Vhr/10Vhr = 0.002 
(a^ + b^)^/^ = 0.0036 
cos a = 3/3.6 = 0.833; a = 33.55° 


From equation (91) 

sin(i(; - \|^^) = Alp = -sin(ip^ + 33.55°) x 3.6x10 ^ 

This dependence of Aip on ipi for -180° < ipj. - 180° is shown in figure 7. 

It shows that Aip can be as large as 0.2°, which exceeds the limitation in 
equation (87). However, the top row in equation (89) can be used to determine 
Alp, that is, 

= [cos(ip - ip^) - cos ^)]^S ^ 

In the example, ip - ip^ = 0.2°. Assuming an even more conservative error, 
ip - ipj_ = 1°; cos(ip - ipi) = 0.99985. By subtracting the known value Qg cos X, 
the residual due to the first term is 

cos(ip - ip.)^^ cos X - cos X = cos X x 0.00015 = -10°/hr x 0.00015 
^ s ^ ^ 

= -0.0015° /hr 
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which is negligible. The remaining terms in equation (92) are, after normal- 
ization by fig cos X 


- fig cos X 

® ^ fi — cos X “ 

s 

Using the same notations as before, we have 

e = -(cos cos ot - sin sin a) (a^ + (94) 

By causing a perturbation on (e.g., a sine wave of known amplitude), 

we have 

6e = (sin cos a + cos sin a) (a^ + b^)^/^6i(;^ (95) 

But from equation (91), it follows that 

6e = (96) 

and 


Alp 


6e 

Sip. 

1 


(97) 


thus determined is added to the value of which nulls fij^ in accor- 
dance with equation (90). The new value + Aij; = is now applied to the 

resolver, which now yields 


1 


. * 


L’r_ 



cos(iJ; - ijj.*) - cos cos X) + sin , *(f|j/$^_ cos X) 

1 Id® X d ® 

-sin(\|j - ip.*) - cos i(;.*(f|j/fi_ cos X) + sin 4'.*(tj/fio cos X) 
X X d ® X d ® 


fig cos X 


(98) 


But 


’I' - = rp - - Alp = 0 


Rearranging and subtracting fig cos X in the first row, we have 


^ * r\ -V t • • I 

- fig cos X = -cos [p^ + sin \p_, n 


i ’'■d 


= -Sin ip^ Cj - cos ip^ 


From the last two equations the drift rates are readily determined 


(99) 

(100) 
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(101) 


"^d" 


— A 

cos 

sin ip.^ 
1 

- ^^g cos A 

_^d_ 


-sin 
L X 

cos ip,* 

A * 

L \ J 


The procedures of gyrocompassing, equations (89) and (97), and calibra- 
tion of equation (101), are indicated in figure 6. 

With the drift ^^te biases thus established, the appropriate torquing 
commands and are applied to the LEG, as indicated in figure 6. 

With the drift rates compensated, the initial gyrocompassing step, defined 
in equations (89) and (90), is repeated with new essentially zero drift rates 
and f)j, so that ^i = ^ is determined with a negligible error and is 
applied as the initial condition i/jq to the integrator in figure 4. The 
computed Earth rate components cos A sin and cos A cos 

are now applied to the torquers T^ and Tg, respectively, as indicated in 
figure 6* Switches and Sy, shown in figure 4, are closed simultaneously 
on application of and This is incorporated in the command mode 

shown in figure 6. 


The azimuth gyro (AZG) is subjected to the Earth rate component Qg sin A, 
Its output is = Qg sin A - ip^,^where is its drift rate. By subtract- 

ing the known quantity, Qg sin^A, is readily determined and the appropriate 
compensatory torquing signal is applied to T^, as indicated in figure 6. 


With the drift rates and eliminated by the foregoing procedure, 
the GIMU is also automatically calibrated, as can be verified by (61) and (62). 
The bias terms g?^ and gq^, which are due to initial offsets, are measurable. 
They can easily be eliminated at the GIMU output, as indicated in figure 6 by 
the calibration signals g?^'^ and gq^^^, respectively. 


Summary of gyrocompassing and calibration procedure: 


1. Open Sjj and Sy and turn on power to LEG and AZG. 

2. ^ After warm-up time, apply and fij^ to ijj resolver, varying 

until f|R = 0. Store 

3. Apply perturbation to and determine gyrocompassing error Aip 

by equation (97). 

4. Determine + Atji. 

5. Apply \p^* to the \p resolver and determine and q^*. 

6. Solve drift rates |(j, q<j by equation (101) . 

7. Repeat gyrocompassing procedure in accordance with step 2. 
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8. Apply compensating torques nc'^, ?c^, r)^^, and d to Tp, T„, 

T^, respectively, simultaneously close and Sy. 

9. Check GIMU outputs Fxp and Fyp. If necessary, null outputs by adding 
&K± 2 nd grii* as indicated in figure 6. 


VI. ERROR ANALYSIS AND EVALUATION 


As shown in section V, errors originate primarily from angular drift rates 
in the LEG and AZG gyroscopes, which, for a given type of gyroscope, are essen- 
tially determined by temperature variations and g-loadlng. A miniature DTR 

gyro, described in appendix E, is chosen as the basis of the error analysis in 
this section. 


Other errors originate in uncompensated disturbances due to friction 
torque, torquer scale factor errors, and gyrocompassing and alignment errors. 

In the analysis, x denotes position deviation, either lateral or longitudinal. 
In the numerical evaluation, the following constants are repeatedly used: 

R = 6.378x10 m, oig = 4.46 rad/hr, = 19.93, 180°/7 t = 57.3. 


A. Sensor Errors 


1. Non-g-sensitive drift 


From appendix E, the average turn-on to turn-on uncertainty is 0.03®/hr. 
This error is calibrated by the procedure described in section V and is com- 
pensated by^the appropriate gyro torquer inputs. Assuming a compensation 
error of 10%, the constant residual drift rate is d = 0.003°/hr. From fig- 
ure 8 in appendix A, for a = 0 at 1 hr U = 4.46), u = 5.4 so that 


0 


X 


^ ^ 6.38x1q6 X 3x 10~3 
cOg 4.46 X 57.3 


5.4 = 403 m 


2. Temperature sensitivity 


Assuming an average temperature sensitivity = 0.03®/hr/®F (appendix E) , 
and the average temperature variation is a constant 0.5°F/hr, the rate of 
change of the drift rate is d = 0.015°/hr2. The resulting positional error 
obeys the same law as given in equation (D-4) . At 1 hr, from figure 10, u = 9 
and the position error is 


^ ^ ^ 6.38x1q6 X 0.015 

m2 19.93 X 57.3 


X 9 = 754 m 


3. g - Sensitivity 

Assuming an overall average sensitivity of Sg = 0.05“/hr/g, known for 
each individual gyroscope (appendix E), and assuming a root mean square (rms) 
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g input, due to turbulence, of 0.25 g, the rms drift rate is 0.0125° /hr. (Due 
to the low-g environment, g^ dependence is negligible.) The bandwidth a 
(defined in appendix A), based on standard turbulence models (ref. 13), is in 
the range of 0.1 to 1.0 rad/sec. Thus, a = a/cOg is 1.5 to 15. An average 
a = 2 is chosen. At 1 hr, u = 5.4 (eq. A-11) 



6.38x1q6 X 0.0125 
4.46 X 57.3 


X 5.4 = 1685 m 


Assuming that Sg is known to 5% or better, this error can be compen- 
sated to l/20th of its value, by applying FjjpS|jj + FypS|y and 

FxpS^x + ^ypSgy to T^, Tp as indicated in figure 6. Thus, the residual 
error is of the order of 80 m and is negligible. 


A sustained acceleration (e.g., as a result of a turn, causing an 
Fyp — 10 m/sec^, for a duration t = 100 sec) can be considered as an acceler- 
ation impulse, since WgT « 1. It excites a drift rate impulse, for example, 
in n 

• 

d = FypS^y = 0.05°/hr/g x Ig = 0.05°/hr 


This disturbance excites a position error (eqs. (82) and (83)) 


Ax 


(gd)T 

2 


(1 - cos (DgX) 


9.81x3600^ X 0.05x100 
19.93 x 57.3 x 3600 cos WgT) 


It reaches its peak value of 150 m at t — Tr/ug — 0.704 hr. Using the pre- 
viously calculated compensation factor, this error can be reduced to ~8 m 
and becomes negligible. 


4. Random drift in LEG gyro 

Assuming an average random drift of d = 0.005'' /hr (appendix E) and a 
correlation time a“l = 0.15 hr, a = a/o)g = 1.5 (appendix A). The value of 
u at 1 hr is, u = 3.1 (fig. 8). The resulting positional error is from 
equation (A-11) 


a 


X 



6.38x1Q6 X 0.05 

4.46 X 57.3 


3.1 = 387 m 


5. Drift of AZG gyro 


, Error in compensating Earth rate- Assuming a constant drift rate 
= 0.1°/hr, for moderate latitudes of (X = 45°), fig cos X = 10°/hr, 

C = cos i|^(0.1°/hr X 10°/hr)/57.3 = cos i/^(0.174°/hr^) . Then, in accordance 
with equation (D-4) and for u = 9 at 1 hr (fig. 10), and, e.g., ijj = 0, 
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Ra 


= 


u = 


6 . 38x106 X 0.0174 
19.93 X 57.3 


X 9 = 874 m 


It follows that the AZG can tolerate a drift rate ~10 times larger than the 
NEG. 


Azimuth transformation error- This error is given by the first term in 
equation (10). Setting = /'J'd dt, the positional error, for example, in 
the north direction, is Ax =//VN4'd dt dt'. For example, if Vn = 200 m/sec 
at 1 hr. 


Ax 


200 X 0.1 X 3600^ ^ 
3600 X 57.3 X 2 


6. Gyro Torquer Scale Factor Uncertainty 

Torquer scale-factor uncertainty generates an equivalent drift rate 

of d = (e^Vj^)/R. For Vjj = 200 m/sec and e,p = 3x10“**, 


3x10“** X 200 X 3600 x 57.3 
3.68x106 


0.00336°/hr 


From figure 8, for a = 0, at 1 hr, u = 5 . 4 . Thus, 


a 


X 



6.38x1q6 X 3.36x10“3 
4.46 X 57.3 


X 5.4 = 453 m 


The effect of scale-factor uncertainty on the AZG is readily shown to be 
negligible. 


B. System Errors 


1. Gyrocompassing Error 

From equation (87) the uncorrected components of Earth rotation are given 
by 


AC = Atf) sin fig cos X 

An = At{» cos ip^ fig cos X 

Assuming a gyrocompassing error (sec. V) of Aip = 0.05®, the rate errors, a| 

and An, at moderate latitudes of fig cos X = 10®/hr, are of the order of 

10°/hr X 0.05°/57.3 = 0.0087®/hr. This is equivalent to a constant drift rate 

d for which a = 0 in figure 8. At 1 hr, u = 5.4, and the positional error 

is by equation (A-11) 
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6.38 x 1 q 6 X 0.0087 x 5.4 
4.46 X 57.3 


= 1172 m 


2. Reaction Torque of the AZG 

The AZG, being subjected to heading rate input ip, exerts a torque Hiji. 

If its spin axis is along X^, this torque is around Y^. For H = 8 g cmsec 
and typical values for ^ of 0.1 rad/sec, this torque amounts to 

= 8 X 0.1 = 0.8 g cm. Since the pendulosity WZ of the GIMU is assumed to 
be 200 g cm, this precession torque is equivalent to an acceleration angle 

e = i)H/(W£) = 0.8/200 = 0.004 rad 

or, an acceleration error of eg = 0.04 m/sec^. Assuming that H is known to 
1% accuracy, this disturbance can be computed, scaled, and subtracted from the 
GIMU output as indicated in figure 6. With these assumptions, the residual 
acceleration disturbance is er = 40 yrad. If, ijj, in a turn persists, for 
example, for t = 100 sec, then, since WgT << 1, this can be considered as a 
velocity impulse exciting a position error (see eqs. (82) and (83)), 

gEj-T 

Ax = sin (jjat . 


It reaches its peak of 


. _ 9.81 X 40x10“^ X 100 X 3600 , 

Ax 7 — — = 31.6 m 

4.46 

at t = ir/2cjg = 0.352 hr. This is a negligible error. 

Assuming that the rms value of ip (due to turbulence) is 0.1 rad/sec, at 
an average angular frequency of a = 0.5 rad/sec, the average value of 
a = a/(i)g is 6.7. The corresponding u in figure 9 at 1 hr is u = 0.9. In 
accordance with equation (C-2) , the positional error at 1 hr is: 

0^ = ROgj.u = 6.38x10^ X 40x10“^ x 0.9 = 230 m 

With the above compensation of 1%, this error is negligible. 

3. Initial Tilt Error 

From equations (82) and (83), the effect of an initial tilt error ni> Ci 
denoted here by results in a positional error. 

Ax = g££(l - cos ajgt)/Wg^ 

Assuming that = 100 prad. Ax at 1 hr is 
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Ax 


9.81x3600^ X lOOxlQ-6 
19.93 


(1 - cos 4.46) = 797 m 


4. Bearing Friction Torque Disturbance 

Gimbal bearing friction is modeled as dry friction exerting a constant 
torque, Tf, due to the pitch and roll motion of the airframe, with respect to 
the gimbal axis frame. The significance of this assumption is that the torque 
signals can be modeled as a random binary process fluctuating between +Tf and 
-Tf, with a zero crossing rate a, determined by the frequency of the random 
pitch and roll motion excited by turbulence. It is known (ref. 14) that 
the autocorrelation function of such a binary process is 

<J)ff(T) = T^2 e-2aM 


The acceleration error induced by this torque noise is e = Tf/W£, so that 


(t) 

ee 



Typical values of Tf are 0.05 g cm; and since W£. is assumed to be 
m = 200 g cm, Oe = 0.05/200 = 250 yrad. 

Assuming again, a = 0.5 sec“^, the value of a is a = la/oig = 13.4. 
In accordance with figure 9, the value of u at 1 hr is u = 0.5, and the 
positional error in accordance with equation (C-2) is: 

= ROgU = 6.38x10^ x 250x10"® x 0.5 = 796 m 


Since the friction model is known, an estimate of e can be provided in real 
time, and a large portion of the friction noise ge can be cancelled by 
subtraction. 

5. Dead Band Due to Bearing Friction 

By the same numerical assumptions as in equation (B-4) , the bearing fric- 
tion causes a dead band of 250 yg in the acceleration measurement correspond- 
ing to 250 yrad. It can be effectively linearized by injecting dither signals 
to the LEG torquers T^ and Tr|. Assuming a sinusoidal dither signal of a fre- 
quency fd = 5 Hz (see sec. V) to provide an angular amplitude of ~10 times 
the dead band, that is, aj = 2.5 mrad, the peak torquing rate r^ is 

r^ = “ 2.5x10"^ x 2Trf^ x 57.3 = 4.5®/sec . 

In accordance with appendix E, this is well within the torquer capacity of the 
gyros considered. In accordance with reference 15, the system is essentially 
linearized. 
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6. Gimbal Torquer Scale-Factor Error 


Error In Schuler tuning- A scale-factor error b^, by in the gimbal 
torquers is mathematically equivalent to the gyro torquer scale-factor errors 
discussed in section VI. In accordance with reference 9 it is assumed that the 
errors in b^, by, denoted et are ~10“‘^. In accordance with section VI this 
would produce an equivalent drift rate of less than 0.001®/hr, and a positional 
error (at 1 hr) of 133 m, which is negligible. 

Proportionality error- The position error is directly proportional to 
^b “ 10“^+. For Vjj = 200 m/sec, at 1 hr, the error is 

Ax = lO"** X 200 X 3600 = 72 m 
which is also a negligible error. 


Summary 

The principal contributing factors to positional error at 1 hr of flight 
are summarized in table II. They are treated as statistically independent, and 
the total error is the resultant rms value. The total rms positional error 
after 1 hr of flight is 2210 m or 1.2 n. mi. 

The individual rms error propagation time histories listed in table II 
are shown in figure 11 up to 2 hr, along with the total rms error denoted by 
It is typical that from these results the initial large tilt error of 
100 prad, assumed in section VI-3, causes a relatively rapid positional diver- 
gence in the first 0.5 hr, whereas the constant drift rate of temperature, 
assumed in section VI-2 causes a relatively rapid divergence after ~50 min. 

An alternative model is to assume a smaller initial tilt error of 25 yrad, and 
a stationary random temperature variation of 1° F and a correlation time of 

^ ^ “ 1 min, for which a = a/tOg = 60/4.46 = 14. The results are shown in 

figure 12. The short-term divergence is considerably reduced, but no signifi- 
cant change is incurred around 1 hr of flight. In a third error model, the 
AZG has a drift rate of 0.05®/hr instead of 0.1°/hr; the initial tilt error is 
as in the first model, and the temperature model is as in the second. The 
results are shown in figure 13. 

It is significant that with all three models the positional divergence at 

1 hr clusters around 2 km. The third model evidently has a smaller divergence 

beyond 1 hr. 


VII. CONCLUSIONS 


It has been shown that the low-g environment assumed in this report per- 
mits a significant simplification in mechanization by dispensing with the 
accelerometers, since moderate size gimbal torquers can provide the required 
torques for specific force measurements. This device constitutes a two-axis 
locally level accelerometer, which is easily integrated with an on-board 
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microcomputer to constitute a complete inertial measurement unit. Rapid gyro- 
compassing is achieved by virtue of the large pendulosity of the gimbal system 
and by strapdown type resolving of the LEG outputs. 

Gyro drift rates are determined as part of the gyrocompassing process. 
Since the essential inertial measurements are derived from one sensor (LEG), 
complex instrument alignment is circumvented. Bias errors due to possible 
correlated gyro and accelerometer noise components are not present. The 
gyros, in particular the LEG, being essentially isolated from airframe motion, 
require only minute torquer activity, so that temperature fluctuations are 
also small. In addition, the low-g environment renders the g and g^ sensi- 
tive errors sufficiently small so that gyros of moderate cost can be used. 

Modem miniature DTR gyroscopes are sufficiently small so that both the 
LEG and AZG can be mounted within the inner gimbal g£. Thus, the complete 
assembly is expected to be comparable in size and weight to a conventional 
vertical gyroscope. The error analysis, based on state-of-the-art sensor data, 
demonstrates that a navigation error of the order of 1 n. mi. at 1 hr of 
flight can be achieved. 
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APPENDIX A 


POSITIONAL ERROR PROPAGATION DUE TO RANDOM GYRO DRIFT RATE 


The specific force measurement error in a Schuler tuned system is given 
in equations (82) and (83). and which are angular deviations from the 
local vertical, consist of constant, random stationary, and divergent terms, 
as indicated in equations (69) and (70). In this appendix, the following uni- 
form notation is used for the drift rates ?d ^d contributing to 

^N- 

d - gyro angular drift rate, defined in geographical coordinates (rad/hr). 
cOg -- Schuler frequency (rad/hr) . 

R - distance of vehicle from the center of the Earth (meters) . 

a” ^ “ correlation time of gyro drift rate process (hours). 

a = a/a)g. 

A 

C 03 g t • 

The drift rate is modeled as a first-order Markov process with the autocorre- 
lation function 




(A-1) 


where = V^^-(d) • The corresponding power density spectrum of d is 


■ “d" “ 

o ^ 


2a 


(A-2) 




The spectrum of positional deviation x as derived from equations (82) 
and (83) is 


d -S2(a2 - s2)((,32 + s2)(w 32 - s2) 

In time domain, Var(x) = is given by 

^O ^o 

where the impulse response hx(t), corresponding to Hjj(S) = 1/S(s2 + 
equation (A-3) , is 


(A-3) 


(A-4) 
0 )g 2 ) in 
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(A-5) 


hx(t) = T (1 - cos OJgt) 


w 


Substituting equations (A-1) and (A-5) in (A-4) and normalizing by 
the double convolution integral is 


t .t 




(1 - cos 0)gT')(l - COS 0JoT")e 


-alx'-T" 


dt' dt' 


“s *^o ”^0 

•'o •'o 


=*rr 


(A-6) 


(1 - COS v')(l - COS v")e ^ ^ dv' dv" 


By the change of variables X = v' - v", equation (A-6) is transformed to 


1=-^ f (1 - cos v’) r 

“s® ^’ = 0 •'v’ 


v'c 


[1 - cos(v’ - X)]e 


-a|x| (_ 


dX)dv' 


- — n 

- I (1 - cos V 
^' = 0 


■>|C, 


aX 


[1 - cos (v' - X) ]e dX 


r 


[1 - cos(v' - X)]e dX 


jdv' 


(A-7) 


The result of the integration of equation (A-7) is 


I = 


— — 7 jaC(3a2 + 2)(a2 + l) + 2(l-e “^) (a^ + 1) [«(a cos ? - sin C) - 1] 

a2(a2 + l)2 I 


+ 2a^ e"“^(a cos C - sin ?) - a^(a^ + 1) 4^sin ? - 
- a^(a^ + cos^ - sin^ 


For 


For 


a 0 , I = (? - sin 5)‘ 


a , 


■K 


3C - 4 sin C + 


sin 2c 


) 


(A-8) 


(A-9) 


(A-10) 


Defining u = (I)^/^, the standard deviation of x, Ojj is 
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Plots of u(a, 5) are shown in figure 8. 



APPENDIX B 


VELOCITY ERROR PROPAGATION DUE TO RANDOM GYRO DRIFT RATE 


Using the same notations and definitions as in appendix A, the spectrum 
of the velocity error v is given by 




2a 


^ (a2 - S2)((0^2 + s2)(to 2 _ s2) 


(B-1) 


In time domain, V^j-Cv) = is given by 

.t -t 


0^2 = J* J* h^(T')h^(x")(|)^jj(T' - T")dT* dx" 


(B-2) 


o "o 


where hy(t), corresponding to H^(S) = 1/ (cjs^ + s2) is 


h^(t) = sin mgt 
s 


(B-3) 


Substituting equations (B-1) and (B-3) into equation (B-2) and normalizing by 
a^^g^, the double convolution integral is 


“s •'o •'o 


? . II -a x'-x" j , j , 

sin o)„x sin o)„x e ' ' dx dx 

o o 


1 , 

2 J j ^ 

tOo^ Jr. Jr. 


Sin V e 


I dv' dv" 


"■s 'o 'o 

By changing variables X = v’ - v", equation (B-4) is transformed to 


(B-4) 


= — ^ f sin v'l f sin(v' - X)e dxldv' 

“s ^’ = 0 (•'v’-? ) 

= — - / ■ f sin v*l r sin(v* - X)e*^^ dX + T sln(v' - X)e dxldv* 
“s ^'=0 y'v'-z •'o J 


(B-5) 
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The result of the integration of equation (B-5) is 


I = 




For 


For 


(a2 + 1)^ 

- 2(a sin ^ + cos 5)e + l| 

a 0 , I = (1 " cos O' 


cos^ Z) 


a 00 




Defining u - the standard deviation of v, 


IS 


= Ro^u 


(B-6) 


(B-7) 


(B-8) 


(B-9) 


Plots of u(Cj a) are shown in Figure 9. 
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APPENDIX C 


POSITIONAL ERROR PROPAGATION DUE TO RANDOM ACCELERATION ERROR 


Using the same notations and definitions as in appendixes A and B and by 
denoting the acceleration error as ge (where 0 is an equivalent tilt angle 
from the local vertical), then the spectrum of the positional deviation x is 


^xx(S) 


2„2 


(a^ 


2a 

S2)(0)g2 + S2)(t0g2 _ s2) 


(C-1) 


where 


Og2 = var(e) 

This has the same structure as equation (B-1) and therefore results in the 
same I as given in equation (B— 6) . Consequently, the standard deviation of 
X, Ojj is 


Ojj = ROgU (C-2) 

A plot of u(^, a) is shown in figure 9. 
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APPENDIX D 


POSITIONAL ERROR PROPAGATION DUE TO CONSTANT HEADING DRIFT RATE ERROR 


From equation (86), the uncertainty in determining the Earth rate compen- 
sating precession command is An = cos X. Assuming perfect gyrocompass ing, 
that is, a precise initial value for ipQ, and a constant unknown azimuth gyro 
drift rate then Aip = ^i<jt, so that the equivalent drift rate due to this 

compensation error is 

C = cos X (D-1) 

The positional divergence due to this error in accordance with equa- 

tions (82) and (83) is 




S3(0)g2 + s2) 


In time domain 


a 


X 



(cos o)gt - 1) + 




(D-2) 


Defining 


u 


4 CJ 


- (1 - cos C) 


(D-3) 


is given by 


a 


X 



A plot of u(c) is shown in figure 10. 


(D-4) 
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APPENDIX E 


TYPICAL PERFORMANCE CHARACTERISTICS OF A DRY-TUNED ROTOR 


Physical data 

- Size 

- Weight 

- Angular momentum 

- Maximum rate input 

- Temperature control 


1 in. diam x 1 in. length 

80 g 

8 g cmsec 
100°/sec 

Integral; sensor + heater 

2. Average drift rate characteristics 

- Repeatability (turn-on to turn-on — non-g-sensitive) 

- g-sensitive drift (including quadrature) 

- Anisoelasticity 

- Torquer scale-factor stability 

- Temperature sensitive drift 

- Torquer scale-factor temperature sensitivity 

- Random drift rate with correlation time of 0.15 hr 

3. Environment operating range 

- Temperature -40®F to 160°F 

- g-capability better than 200 g 

- Maximum transient 

rate input 300®/sec 


GYROSCOPE 


0.025°/hr la 

0.02° /hr la 

0.02°/hr/g2 

la 

50^200 ppm 

0.03°/hr/°F 

la 

0.035°/o/°F 

la 

0.005° /hr la 
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TABLE I.- SYSTEM CONSIDERATIONS IN GIMBALLED AND STRAPDOWN MECHANIZATIONS 

SUBJECT TO ASSUMPTIONS 1 THROUGH 9 


Category 

Meehan 

ization 

Gimballed 

Strapdown 

Specific force sensors 

At least^ two in the 10*^ 
dynamic range class 

At least^ three in the 
10^ dynamic range class 

Gyroscopes - dry tuned 
rotor (DTR) - TDF type^ 

At least two^: One in 

the 0.1® /hr class and 
one in the 0.01° /hr 
class. Dynamic range: 
1Q3. 

At least^ two: Both in 

the 0.01®/hr class. 
Dynamic range: 40x10^. 

Attainability of 0.01“/hr 
gyro performance 

Moderate cost due to low 
gyro torquer activity. 

High cost due to high 
gyro torquer activity. 

Mechanical assembly 

At least three gimbals 
with three torquer - 
pickoff pairs. 

Fixed mounting. (Three 
dithering torquers for 
laser gyros.) 

Electronic circuitry 

Sensor interfaces, plat- 
form servos, temperature 
controllers. 

Sensor interfaces, tem- 
perature controllers, 
(control circuitry for 
laser gyros). 

Computational load 

Integration of naviga- 
tion equation, compen- 
sating torques for Earth 
rotation. 

Integration of naviagtion 
equation, compensating 
torques for earth rota- 
tion, direction cosine 
matrix, transformation to 
geographical coordinates. 

— » — 1 

Susceptibility to vibra- 
tion and shock 

Low, due to rotational 
isolation by platform. 

Potentially high, due to 
possible correlated noise 
in gyros and specific 
force sensors. 

Alignment and bias 
calibration 

Insensitive to airframe 
orientation. One cali- 
bration valid for arbi- 
trary flight profile. 

No special devices 
required. 

Sensitive to airplane 
orientation. Many cali- 
bration points are 
required for arbitrary 
flight profile. Special 
calibrating device 
required. 

Availability of attitude 
rate information 

By differentiation of 
gimbal pickoff readings. 

Directly available from 
the body mounted gyros. 


^"At least" refers to possible redundancy requirements. 

^This type is currently a valid comparison with strapdown mechanizations. 
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TABLE II.- RMS POSITIONAL ERRORS AFTER 1 HR 


Section 

Description of error source 

rms error at 
1 hr, m 

1 

VI-Al 

Non-g-sensitive gyro drift 

403 

2 

VI-A2 

Drift rate, due to temperature variation 

754 

3 

VI-A4-3 

Drift rate, due to error in Earth rate compensation 

874 

4 

VI-A4-b 

Distance-related error due to AZG drift 

628 

5 

VI-A5 

Random drift rate in LEG 

387 

6 

VI-A6 

Gyro torquer scale-factor uncertainty 

453 

7 

VI-Bl 

Gyrocompass ing error 

1172 

8 

VI-B3 

Initial tilt error 

797 

9 

VI-B4 

Gimbal bearing friction torque 

796 



Total rms error 

2210 
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Figure 1.- Description and definition of the axis system of the two-gimbal 
gyro pendulum and principal loop closures. 
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Figure 6. Integration of the GIMU in a complete inertial measurement system. 
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Figure 8.- Normalized positional error propagation as a function of normalized 
time c = o)gt with gyro noise bandwidth a = awg as a parameter. 
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Figure 9.- Normalized velocity error propagation as a function of normalized 
time ? = ojgt with gyro noise bandwidth a = awg as a parameter, or 
normalized position error with accelerometer noise bandwidth a = ato with 
a as a parameter. ® 
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Figure 10.- Normalized error propagation as, a function of normalized time 

t, = (Ogt for a constant AZG drift rate or a constant rate of change of 

temperature of the LEG. 
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Figure 11.- Positional 
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error propagation with a 100 
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Figure 12.- Positional error propagati 9 n with a 25 prad initial tilt error and 
constant drift rate of the AZG of = 0.1“/hr and a stationary random 

temperature variation of the LEG of 1° F. 
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Figure 13.- Positional error propagation with a 100 yrad initial tilt error a 
drift rate of the AZG of = 0.05 “/hr, and a stationary random tempera- 
ture variation of the LEG of 1® F. 
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